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Abstract 


The flux of energetic protons In Saturn's inner magnetosphere was 
observed in two channels from 48-63 and 63-160 MeV. Absorption features due 
to the G ring and the satellites Enceladus and Mimas were easily identifiable. 
The flux observed in the absorption slot of Mimas can be maintained by the 
decay of a cosmic ray albedo neutron flux of 7 x 10“^ cm’^ s"^ sr’^. This 
flux is entirely consistent with calculations of the neutron flux produced by 
galactic cosmic ray interactions with the rings of Saturn. The omni- 
directional oroton flux of 8.2 x 10^ cm"^ s“^ at 2.734 requires a residence 
time of 30 years. Both the residence time and the energy spectrum are 
comparable to those found in the inner radiation belt of the Earth. The 
angular distribution is nearly isotropic in the Mimas slot and beyond 4Rg. 
Otherwise the pitch angle distribution is pancake and has been approximated by 
sin^e with n in the range 2-7. This distributior is consistent with an 
isotropic neutron source in the ring plane. 
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Introduction 

Intense penetrating radiation was discovered by Pioneer 11 In Saturn's 
Inner magnetosphere. This radiation was found between the orbit of Mimas 
(3.075 Rg) and the outer edge of the A ring (2.265 R 5 ). Protons with energies 
above 80 MeV are a major component of this radiation (FIlHus et al., 1980; 
Fniius and Mcllwaln, 1980; Van Allen et al . , 1980). Initially, we did not 
fully appreciate the sensitivity of the Goddard Space Flight Center/University 
of New Hampshire low energy detectors to the Intense flux above 50 MeV 
(McDonald et al., 1980; Tralnor et al., 1980) and attributed our Pioneer 11 
obssrvat’ cns to 1.1 - 8.1 MeV protons. Subsequent more detailed analysis of 
the response of these detectors showed that they could have responded to 
either low or high energy protons but that based on the observed flux above 
50 MeV the contribution of high energy protons would be substantial inside the 
orbit of Mimas. A similar conclusion was reached by Simpson et al. (1981). 

Mimas should stop the slow Inward diffusion of protons, but McKIbben and 
Simpson (1980) showed that the rapid episodic diffusion past Mimas during 
interplanetary disturbances was consistent with their initial interpretation 
of the data. Another explanation is based on the decay of energetic neutrons 
produced by cosmic ray Interaction with Saturn's atmosphere and rings. This 
GRAND source is responsible for the high energy protons In the Inner radiation 
belt of the Earth, and estimates indicated that a trapping lifetime of ~ 20 
years would be required at Saturn to permit the flux to build up to the 
observed values (Fllllus and Me II wain, 1980; Van Allen et al., 1980). 

Voyager 2 carried two experiments Identified as the Cosmic Ray Subsystem 
(CRS) and the Low Energy Charged Particle experiment (LECP), which observed 
the protons In the inner magnetosphere. Preliminary Interpretation of the new 
data showed that any change In the high energy proton flux between 1979 and 
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1981 fen within the uncertainty of the measurements (Krimigis and Armstrong, 
1982; Vogt et a1., 1982). The LECP experiment also reported the proton 
spectrum and found a low energy component (E < 0.5 MeV) In addition to the 
high energy flux (Krimigis and Armstrong, 1982). The three-energy channel 
LECP measurement of the high energy proton flux peaks In the 54-87 MeV channel 
with lower fluxes near 10 MeV and 100 MeV. Good agreement was found In the 
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absolute flux above 63 MeV between the CRS and LECP measurements. 

The CRS on Voyager 2 was able to measure the proton flux in the Inner 
magnetosphere In two channels (48-63 MeV and 63 - 160 MeV). The counter 
telescope used for this measurement a low electron detection efficiency 
and was sensitive enough to collect adequate statistics In the Important 
absorption slot produced by Mimas. This constitutes a major advantage over 
other experiments; however, the greater sensitivity also Implied substantial 
dead- time and accidental coincidence corrections. 

The Voyager 2 trajectory In the Inner magnetosphere covered the latitude 
range from 30“ to -20“ . The latitude difference between the Inbound and 
outbound passes permitted a measurement of the flux at two points on a field 
line, and Intensity differences can be related by Llouvllle's theorem to the 
proton flux at two different equatorial pitch angles. These mea'surements do 
not permit an Identification of the shape of the angular distribution; still 
they add to the Pioneer results which were either obtained with omni- 
directional detectors or shielded detectors whose angular response was 
uncertain because of shield penetration. 

The purpose of this paper Is to present a more careful analysis of our 
observations than was possible In the preliminary publication (Vogt et al., 
1982) and to Interpret the observed proton fluxes in terms of a balance 
between the proton source and the absorption due to the G and E-rIng and the 
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satellites Enceladus end Mimas. We derive a new expression for satellite 
absorption for the case when the gyroradlus of the particle Is larger than the 
radius of the satellite and explore the effects of absorption by a satellite 
or dust ring on the energy spectrum of Interacting protons. The analysis 
permits us to deduce the flux of CRANO neutrons required to maintain the 
observed proton flux. The strength of the neutron source Is consistent with a 
source at the rings but too high for an atmospheric source. Our results do 
not find a natural explanation In terms of episodic diffusion. 

Instrumentatl o n 

The CRS experiment has been described by Stone et al. (1977), and the 
data used in this paper were collected with the High Energy Telescope (HET). 
The high energy proton mode of the HET has three proton energy channels. 
Protons enter through detector B1 (Fig. 1), and a coincidence-anticoincidence 
matrix determines whether or not they stop In detector C4, C3, or C2. Alpha 
particles and heavy Ions are discriminated against by an anticoincidence with 
a slant threshold. Electron sensitivity Is minimized by high thresholds on 
detectors B2, C4, C3, and C2, which are well above the average energy 
deposited by an electron. Thus the higher order coincidences have very low 
electron sensitivity even though complete proton-electron separation Is not 
possible In an Individual detector. 

The primary function of the HET Is to monitor cosmic rays; consequently, 
the electronic circuits were designed with long time constants for good 
stability and low power consumption. At high rates, the anticoincidence 
circuits block only a fraction of the pulses and accldentlal coincidences 
become a problem. In order to operate In Saturn's Inner magnetosphere, we 
turned off the low threshold detectors Bl, Cl, and the guard counters (Fig. 1, 
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Table 1). This offers two advantages: first, the efficiency of the 
anticoincidence shield becomes Irrelevant and secondly, the solid angle for 
true coincidences Is Increased from 0.96 cm^sr to 10.5 cm^sr for B2, C4, C3. 
and C2 coincidences. Because the accldentlal coincidence rate remains 
unchanged, the true- to accldentlal -coincidence ratio becomes much more 
favorable. With this change the solid angles and proton energy ranges are 
less well defined and corrections are required for particles entering or 
escaping through the sides of the telescope. In addition, most alpha 
particles are now counted In the proton channels. They can no longer trigger 
the slant threshold because the B1 contribution to the summed pulse heights Is 
absent. The sensitivity to alphas should present no problem In the Inner 
magnetosphere where protons dominate the energetic particle population 
(McDonald et al., 1980; Simpson et al., 1980), but does change the cosmic ray 
response. 

Corrections to the raw data were performed In two steps. The true 
coincidence rates were determined by correcting for accidental coincidences 
and losses due to the dead time of the electronic circuits. In the second 
step, we calculated the flux In specific energy ranges by taking account of 
the geometric factor and correcting for higher energy protons going through 
the sides. These corrections depend on the shape of the proton spectrum and 
were carried out relteratively for power law and exponential spectra until 
neither the spectral Index nor the flux changed. Errors In the final results 
were calculated by introducing appropriate uncertainties Into the different 
corrections and adding all errors in quadrature. The flux In the 27-48 MeV 
range can. In principle, be derlwtd from the B2, C4, CT rate; however, the 
corrections for higher energy particles passing through the side were so large 
that no reliable flux could be obtained In this energy range. A further 
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description of the logics and how the corrections were calculated Is given In 
the Appendix. 

Observatl ons 

The corrected coincidences rates are shown In Fig. 2, where R£ represents 
the B2, C4, and C3 coincidence rate with C2 In anticoincidence and R 3 
represents the B2, C4, C3, and C2 coincidence rate. Since the two rates are 
almost equal, the R 2 rate has been shifted down by 1 decade. Both the R 2 and 
R 3 rates were nearly constant from 10 to 5 R 5 . They Increased by three orders 
of magnitude and reached their maximum shortly after perlapsis when Voyager 
was near the tnagnetic equator. Strong absorption features are seen at the 
orbits of Enceladus and Mimas. During the Inbound pass, a spacecraft slew was 
performed near the G ring and produced a dip In the counting rates when the 
detector was pointed parallel to the magnetic field direction. The absorption 
of the G-rIng at a constant pitch angle was observed during the outbound pass 
and produced only a change in the slope of the proton Intensities. 

The HET detector was reconfigured to the special operating mode at 
10 Rg. Between there and 5 Rg, Wie observed coincidence counts were primarily 
due to cosmic ray protons and alpha particles. This deep penetration would 
not be expected on the basis of the Stormer cut-off rigidity (Table 2); 
however, both Pioneer 11 and Voyager 1 observed that the actual cut off of low 
energy protons (2-10 MeV) occurs between 7.5 and 10 R 5 (McDonald et al., 1980; 
Vogt et al., 1981). It Is believed that these low-rigidity particles gain 
access to the magnetosphere via the magneto tall. The presence of alpha 
particles can be deduced frcrni the summed pulse heights of the C2, C3, C4 
detectors because an appreciable fraction of the events has pulse heights In 
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excess of the total energy that can be deposited by a proton. Thus the 
trapping boundary for > 48 MeV protons may not be too far beyond 5 Rg, 

The intensity increase starting inside the orbit of Tethys can be 
attributed to energetic protons. A more detailed examination is required to 
identify the source of the counts observed in the Mimas absorption feature 
between 3.02 and 3,14 Rg. A careful examination of our corrections showed 
that these counts cannot be due to accidental coincidences. The minimum 
dipole cut-off rigidities in this region are so high (Table 2) that even a 
somewhat reduced cut-off due to magnetotail access would filter out the large 
majority of the cosmic rays. The R 2 and R 3 counting rates are about the same 
at 3,08 K$ as between 5 and 6 therefore, residual cosmic rays can account 
for only a small fraction of the observed rates. We are either observing true 
proton coincidences or a residual electron sensitivity. Whether or not 
electron coincidences are present can be determined on the basis of the ratio 
between the singles rate in detector C3 and the particle flux derived from the 
coincidence rates. The geometry is such (Fig, 1) that detector C3 is shielded 
relatively uniformly with a threshold of approximately 48 MeV. If the coin- 
cidence counts are due to protons, then the flux above 48 MeV derived from the 
coincidence rates should be approximately equal to the C3 rate divided by its 
omni-directional geometric factor of 56.6 cm^sr. This was true within ± 30 
percent in the Mimas gap as well as at other places where we were clearly 
observing trapped protons. In contrast, if the coincidences were due to 
electrons, then the C3 rate should give a flux ^ 20 times larger than the 
coincidence rate because the low detection efficienty for electrons in 
detectors B2, C4, and C2 would give a much lower coincidence rate for equal 
electron and proton fluxes. This demonstrates that the coincidences are 
indeed due to a residual proton flux in the Mimas gap. Our data are 
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consistent with no electrons reaching C3; however, we cannot rule out the 
possibility that - 1/3 of the singles counts in C3 are due to electrons. 

Fig. 3 gives the proton fluxes between 48-63 MeV and 63-160 MeV corre- 
sponding to the rates shown in Fig, ?., The flux was derived from the 
corrected rates by dividing by the geometric factor and correcting for 
particles going through the sides (see the appendix). The corrections in the 
48-63 MeV channel were much larger than for the higher energy channel and 
depended on the slope of the spectrum. The flux in the two channels and the 
slope were calculated reiteratively. The values shown in Fig. 3 were 
calculated for an exponential energy spectrum (j « exp - E/Eq^), but essentially the 
same values were obtained for a power law spectrum (j « E“T), Except during 
the spacecraft slew near the R-ring, the telescope pointed almost perpendicu- 
larly to the local magnetic field (Fig. 4). Near the Mimas absorption region, 
the 48-63 MeV flux could not be recovered because after correction the errors 
were comparable to the flux. The error bars in the inner magnetosphere are 
not due to counting statistics but reflect the uncertainties of the various 
corrections; thus, the curves are smooth in spite of the large error bars. 

The spectral indices derived from the fluxes of Fig. 3 are shown in Fig. 

5. No change in either Ech T is noticeable at the orbit of Tethys; however, 
a distinct softening may be noticed near the orbit of Enceladus. We have no 
data on the spectral parameters in the Mimas absorption gap. The spectrum 
hardens at the G~ring (2.8 R 5 ) and appears to soften again at 2.734 Rs. 

Although the absolute errors are large near periapsis, relative values should 
be significant. The value of Eq^ - 40 MeV observed at 2.734 R$ is similar to 
spectra found in the earth's radiation belt at L 2.1 (Lavine and Vette, 

1970). 
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In the inner magnetosphere identical fluxes at equal distances were 
observed between the inbound and outbound passes by Pioneer 11 but not by 
Voyag-'r 2 (Figs. 2 and 3). The difference is that Pioneer stayed very near 
the equator while the latitude of Voyager 2 changed from +30® inbound to -20® 
outbound. The flux ratio at equal Equatorial distar<*es between the two 
Voyager 2 passes provides information about the pUch-angle distribution. The 
calculated Equatorial distances depend critically on the magnetic field model 
used. As shown by Vogt et al. (1981, 1982), a centered dipole model is 
inadequate because the satellite absorption features are off-set relative to 
the position of the satellite. In this analysis we used the '73" model 

^uuiiMci Mcy cu □!., iyoc) uu uDuaiii uiic iieiu &(.rengbii onu cqUdburidi Crossing 

of magnetic field lines going through the spacecraft. Tables of this model 
were kindly made available by Drs. Ness, Acuha, and Connerney. As can be seen 
from Figs. 2, 3 and 7, this model places the absorption features at the 
expected positions. 

The experiment measured the integral flux J(R, e^) at a local pitch angle 
0L and field strength Bl- The local field and pitch angle are related by 
magnetic moment conservation (Schulz and Lanzerotti, 1974) to the Equatorial 
pitch angle 0 and field strength B by: 


sin^e = (B/Bj^) sin^e^ ( 1 ) 

Equatorial pitch angles corresponding to 0j_ are shown in Fig. 4. According to 
Liouville's theorem, the local and Equatorial flux are the same at related 
pitch angles. We approximated the angular distribution by sin'^0 to obtain the 
Equatorial flux at 0 = 90° as a function of the local flux. 
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3(R, 9Q0) = 


J(R, 6|_) 


(2) 


where U(R, e) is the Equatorial flux and R is the radius at which the field 
line crosses the Equator. 

The use of sin'^0 for the angular distribution is an ad hoc assumption. 

The non-gyrotropic or first order anisotropy, which is not included in this 
expression, should be negligible because it depends on the ratio of the 
corotation velocity to the particle velocity and on the relative change of the 
flux over a distance of two gyroradii. Both of these effects are small except 
possibly for a gradient effect at the edges of the Mimas gap. A better 
approximation to the pancake pitch-angle distributions in the inner 
magnetosphere would be "a + sinf'e" for 0 in the range of 0 to 180°; this 
avoids the singularity at small pitch angles. Cecause we can derive only one 
parameter, we set "a" equal to zero to obtain a lower limit on the value of 
■'n". The sin^e angular distribution is not expected to hold in satellite 
absorption regions because particles are absorbed preferentially near 90° and 
a cos 0 • sinf^0 distribution would be more appropriate. A negative value for 
n indicates that this is happening. 

The value of n can be calculated from the intensity ratio between the 
inbound and outbound passes and is given by: 


log (J^n'^'^out^ 
Tog (sin 6,^/sin 


( 3 ) 


where and 6out ai"® the Equatorial pitch angles. It should be noted that 
the effective angular resolution of the detector is improved berauco a nix . — 
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The result shown in Fig, 6 does not represent an actual angular distribution 
but does indicate how the distribution changes with distance. The 

measurements are consistent with an isotropic distribution outside of 4 R$» 
but the observationc are relatively insensitive because they were made for 
equatorial pitch angles below 45®. Inside of 4 Rs» n increased rapidly to a 
value of 7 at 3.26 Rg, This is the most anisotropic distribution we 
observed and occurs near the -inner edge of the E-ring (Baum et al,, 1981). 

The large anisotropy indicates that ring absorption must have been small "^n 
that region and is consistent with the inability of Voyager 1 to image the E 
ring inside of 3.5 Rg (Smith et al. 1981). 

On either side of the Mimas gap, angular distributions may be approxi- 
mated by sin^e. In the absorption gap, our measurements indicate a dumbbell 
distribution, but an isotropic distribution still falls within our errors. 

The spacecraft slew at - 2.8 Rg gives an independent observation of the pitch 
angle distribution of the proton flux in the G ring which can be represented 
j/y 0.4 + sin e. The actual angular distribution must be more anisotropic 
because of the loss of resolution by the - 70® acceptance cone of our 
detector. Inside the G ring, the distributions become more pancake again, but 
the closeness to the equator and large acceptance cone of the detector 
preclude meaningful measurements. In this region. Van Allen et al, (1980) 
found a sin^e distribution and Krimigis and Armstrong (1982) reported sinBe.^ 

A phase space density analysis was performed to identify the proton source 
region and to determine the direction of diffusive flow from the slope of the 
density. For non-relativistic particles, the Vlasov distribution function in 
momentum space is related to the differential flux j by f = j/(2m E), where m 
is the proton mass and E its energy. Because our measurements give only the 
integral flux, J, between lower, Ei, and upper, E£» energy thresholds, we 
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have to use the shape of the spectrum to calculate the distribution 
function. We also extrapolated our measurements to an equatorial pitch angle 
of 90° in order to follow the change in density of protons that mirror near 
the equator. The following expressions were used: 


f{E, R. 900) 


J(R. 6l) 


g(E) 


(4) 


g(E) 


(y-1) 

‘'i 


E"(y+1) 

" h 


(5a) 


exp - E/Eqi^ 

' E Ej„t(exp - Ej/Ejh) - (exp - 

where equ. (5a) applies for power law spectra (j « E"T) and eq. (5b) for 
exponential spectra (j « exp - E/Ech)* Values of the spectral index y or 
characteristic energy Eq^ gre shown in Fig. 5, 

We want to follow particles with a constant first invariant that is with 
a fixed magnetic moment as they diffuse inward; therefore, the distribution 
function for non-relativistic particles has to be calculated at energies 
related by: 


E= 42 -B(R) (6) 

0 

where Eq and Bq are the reference energy and magnetic field and B(R) gives the 
equatorial field at distance R, The magnetic moment was chosen such that the 
proton energy fell primarily into the 48-160 MeV range. This restricted the 
analysis to 15 BeV/6 particles inside the orbit of Enceladus and 60 BeV/G 
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between 6 and 4 Rs- 7 shows the results for an exponential spectrum and 

characteristic energies of 40 MeV for the 15 BeV/G curve and 33 MeV for the 60 
BeV/G curve. 

As can be seen from equ. (4), the value of the distribution function 
depends on the angular distribution through the parameter n of the sin^e 
distribution. In order to demonstrate that the angular distribution changes 
with distances, we plotted the relative phase space densities based on a sin2e 
distribution in Fig, 7a; Fig. 7b shows relative phase space densities for the 
angular distribution parameters shown in Fig. 6. A comparison between Figs. 

7a and 7b shows that the major trend versus distance of the distribution 
function is not changed by specific assumptions about the angular 
distribution. 

Relative phase space densities calculated with power law spectra are very 
similar to those shown in Fig. 7. A minor difference occurs between 4.2 and 
5 Rs for 60 BeV/G protons; the power law spectrum with y = 2.2 gives a 
noticeable minimum at 4. 7-4.9 Rs« The density ratio between 3.3 and 4.3 R$ 
for 15 BeV/G protons is 30 for = 40 MeV but only 11 for y = 1.5. Such 

differences are to be expected because the radial dependence of the 
distribution function is sensitive to the slope of the energy spectrum. For 
the purpose of this paper, however, these Inaccuracies are not important 
because they do not alter the conclusions to be discussed later. 

The slope of the distribution function shown in Fig. 7b vanishes at 2.8 - 
2.90, 3.33, and in a broad region between 4.5 and 5 R 5 . A slight dip or at 
least a flat region occurs at - 2.8 Rs» where both the first and second 
derivatives vanish. On this basis, the diffusive flow of protons between 3.33 
and 4.5 R$ is towards Enceladus from both sides, and the flow between 2.90 and 
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3.33 Is towards Mimas, Protons from inside 2.8 R$ diffuse on the average out 
towards the 6-ring, but the the diffusive flow vanishes in the G ring. 

Proton Absorption by Moons and Dust Rings 

The absorption of trapped particles has been treated by several 
authors. We follow the approach of Mead and Hess (1973) and Thomsen et al. 
(1977). A characteristic lifetime against absorption is calculated by 
averaging over gyroperiod, bounce period, drift period and the range of radial 
distances between the absorber and the gyrocenter of the particle. The 
geometric absorption cross section is used for the satellite. This approach 
does not pick up minor resonances such as between drift period and bounce 
period; however, these resonances are in general too narrow to contribute 
significantly to the lifetime. Previous work was primarily applicable to 
lower energy particles. The effective absorption cross section of the 
satellite for these particles is that of a sphere with a radius equal to the 
satellite radius plus the gyroradius of the particle. All particles are 
absorbed if the motion of the gyrocenter intercepts the sphere. This simple 
model underestimates the characteristic lifetime of high energy particles if 
they are absorbed by a small satellite like Enceladus or Mimas. Since our 
initial writing of this paper, we have become aware of an internal University 
of Iowa report (Rairden, 1980) which treats this problem numerically. Their 
discussions and numerical simulation support the analytical treatment given 
bel ow. 

Let us first calculate the characteristic lifetime of protons which are 
injected uniformly over the sweeping region of a satellite. This is the case 
for a GRAND source and the eccentric orbit of Mimas. Protons drift relative 
to the satellite with a period xd» ^nd the fraction P is absorbed each time they 
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drift past the satellite. An average value of P is required because P is a 
function of the distance between Wie satellite and the path of the gyrocenter. 
If j is the differential proton flux at energy E, then <P> • j Is the decrease 
in flux per drift period. For a stable proton population, we need a source Sp 
of such strength that: 


SpT^ = <P> • j (7) 

If the proton flux is given in units of cm“^ s"^ sr"^, the Sp is in units of 
cid"2 s” 2 sr“^. This permits us to define a characteristic proton lifetime 
against absorption 


""p " ^ 


( 8 ) 


such that it gives Sp = j/tp. In the following, we use the equations of 
Thomsen and Van Allen (1980) to describe the particle motion but use the 
equatorial magnetic field at Mimas of 0.007338 gauss given by the Z3 model. 

The equation for the particle bounce period contains the integral H(e) 
along the field line, and as Birmingham (1982) has shown, the dipole value of 
this integral is still applicable at 3 Rg. 

M 0 i*t of the symbols to be used are defined in Table 3. The lifetime of a 
proton population on a specific drift shell depends on the probability that 
protons can drift past the satellite without being absorbed. This can happen 
because a) a proton may be near a mirror point and does not cross the equator 
as it drifts past the satellite; b) even though It crosses the equator it may 
be in the wrong gyrophase to hit the satellite; and c) for an eccentric 
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satelUte orbit, the satellite may be temporarily too far away from the drift 
shell to be Intercepted. 

Each particle crosses the equator twice during a bounce period; 

therefore, 2tg^/-tjj gives the fraction of the particles that cross the equator 

during the time tg to be derived shortly when their trajectory Intercepts the 

satellite. If the gyroperiod Xg < tg, a particle with a 90** pitch angle 

cannot drift past the satellite without being lost; however. If © 1= 90®, It 

may skip the satellite vertically provided It moves more than the satellite's 

diameter in one gyroperiod (Fig. 8c). For a particle velocity V, the vertical 

distance is given by V t , and the fraction lost is 2r/V.T_, provided that 

By ■ 9 

r/{TTp) « 1. As McDonald et al. (1980) have shown, vertical skipping reduces 
absorption probabilities considerably for > 2 MeV protons at Mimas. In an 
eliptic orbit, a satellite spends a larger fraction of Its time at either 
extreme than near the middle of the absorption region; however, for the 
purpose of this calculation we have assumed equal probability for any distance 
between the maximum and minimum radius of the orbit. In that case 2{p + r)/AR 
gives the probability that the gyrocenter Is close enough for absorption. 
Combining the three factors, we find the probability that the particle Is 
absorbed as Its gyrocenter drifts past the satellite: 


(„ 1 ) 


2r 

tgcos 


^ /2{p + r)^ 

'6 ' aR — ^ 


(9) 


In equation 9, V. has been replaced by V cos©. It should be noted that the 

n 

absorption probability due to any one of the three factors discussed above 
cannot exceed unity, and expressions given In the three brackets of equ. 9 
become Incorrect when this happens. In particular, the second factor should 
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be set equal to I'for cosB < 2r/(VTg)» that is near 90®; the third factor is 1 
at Enceladus where AR < 2(p + r). 

The period of time during which particle trajectories intercept the 
satellite can be calculated from the encounter geometry. As shown in Fig. 8, 
the gyrocenter has to be located somewhere in the interval of length D if the 
leading side of the gyrocircle is to intercept the satellite. An equal 
distance is available for an intercept by the lagging side. If v^j is the 
drift velocity, then tg is given by 


t 

a 


2D 0 Tj 
v^ = — 


(10) 


v^ has been expressed in terms of id in the second part of Equ. 10. D is 
given by the following expressions (Figure 8a, b): 


D = /(p + r)^ - x^ - /(p - r)^ - x^ for x < (p - r) (11) 

= /(p + r)^ - x^ for (p - r) < x < (p + r) 


where x stands for the distance between the center of the satellite and the 
drift path. The average value of D with equal weight for each x is 

<0> = C/o O'!*] / iV 

<D> = . T . for p > r (12) 

p + r 

= ^ (p + r) for r > p 


and 


ORIGINAL PAGE IS 
OF POOR QUALITY 


17 


P » 

" FTF^TT 


for p > r 


(13) 


If we substitute Equation (13) for <t^> and » 4R H(9)/V In Equation (9)* we 
find the average absorption probability per drift period: 


<P> 


2 sine 

AR Tg H(0) cose 


(14) 


H(e) is the line integral between mirror points, defined in Equation 14 of 
Thomsen and Van Allen (1980) and assumes values between 0,74 at e « 90“ and 
1.38 at 0 * 0“; pQ is the gyroradius for 8 ® 90“ . 

From Equation (8), (14). tg/PQ * 2it/V and the relation between particle 
velocity and energy, we find the characteristic lifetime at equatorial pitch 
angle e: 


AR H(8) . H(6) . J JJ o 

P V r^ tan 0 ,c r^ tan 0 /t(h, + ’dt^) 

for 0 * 90® or 0® 

T = 6.988 X 10^ ^ — for protons at Mimas 

p 

The particle energy E and rest mass Eq MeV, c is the velocity of light. 

and Tp is in seconds. Equation 15 does not hold at 90“ where it 

underestimates tp, nor near 0“ (pq sin 0 < r) where it overestimates xp« 

However, we can find the lifetime averaged over pitch angle <x >. from 

P 0 

Equation 15 if the solid angle subtended by these regions is small. For an 
isotropic pitch angle distribution, this gives at Mimas: 
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<T “ 6.99 X 10" 
P 0 


E + 938.21 



(16a) 


Equ. 15 and 16a give relativistic expressions; however, we win use primarily 
the non-relativistic approximation In our energy range E « 938 MeV. In 
contrast to the lower energy expressions {Thomsen et al., 1977), xp Is 
Independent of the drift period, depends sensitively on the equatorial pitch 
angle, and has a weak energy dependence, proportional to E“^^^ for non- 
relativistic particles. 

The above expressions apply to differential energy spectra. To handle an 
Integral flux, we have to average equ. 15 over energy to find <Ve 
alternately find the effective energy at which tp is equal to the average 
value. If J Is the Integral -Flux between energies Ej and E 2 , then 


J 



T S dE “ '^T ^r- 
p p p F 


/e! Sp'lE 


(16b) 


Now we can make use of the non-relativistic approximation to equ. 15 to 
calculate the effective energy Egff which 1 s given by: 

^eff t .'eJ / k\ 

For this calculation ite approximated the source spectrum by the proton 
spectrum at 2.73 R 5 which Is of the forum Sp « exp-E/EQ|^. For the range 63- 
160 MeV, is equal to 90 MeV if Eq|^ ~ 40 MeV; Eg^^ does not depend 

sensitively on the value of Eq^. In addition, xp does not vary rapidly with 
energy; thus, these approximations have only a secondary effect on <xp>£. 

The approach used for uniformly Injected protons has to be modified for 
protons diffusing Into the Mimas absorption region because the gyrocenters are 
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not unifonnly distributed over x but start out at the edges of the absorption 
region. If diffusion is slow, protons are absorbed only every 22 h 37 m when 
Mimas is near the edge. Most of the absorption should occur where the absorp- 
tion probability is highest, that is at x = p - r where 

D = 2 rw (17) 


Using this value for D, we can calculate the absorption lifetime while Mimas 
is at the extremes of its orbit. In that case, only the first and second 
factors in equ. 9 are relevant: 


n Tg H(e) cose 
2r (p^r sine) 


( 18 ) 


T = 8.801 X 10^ 


F f 938.21 H(e) cose 

[E{E + 1876.42)]^^^ /sTnF 


for protons at Mimas 
(E in MeV) 


This lifetime is given as the last entry in Table 3. The 48 min lifetime of 
90 MeV protons at 6 = 85° shows that near equatorial protons are indeed lost 
each time around because Mimas spends each period - 3 h near apoapsis and near 
peri apsis. The distance most of these protons can diffuse in one period of 
Mimas is expected to be small compared to the width of the gap; thus, these 
protons cannot penetrate into the gap. In contrast, at 8 = 45° the proton 
lifetime is 8.8 hrs and about 70 percent survive each encounter; therefore, 
they can penetrate some distance. Thus their lifetime may actually be closer 
to the lifetime of randomly injected protons. These conclusions are 
.consistent with our data. Although the counting statistics are marginal. 
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there is a distinct indicate' on (Fig. 2 and 3) of a square Mimas absorption 
slot for 0 = 75® (outbound) and a rounded slot for 6 = 45® (inbound). 

The characteristic lifetime at Enceladus is also given by equ. 18. The 
orbital eccentricity of Enceladus is small, and the third factor of equ. 9 is 
unity when AR < 2(p + r). Because the width of the sweeping region is less 
than a gyroradius, most of the protons are diffusing towards Enceladus, rather 
than being injected into the gap; thus, equ. 17 gives the appropriate value 
for D. The energy dependence of t is proportional to E-1/^. Consequently, 
absorption by Enceladus will produce a somewhat softer spectrum than the 
injection spectrum. 

Absorption by a dust ring affects the trapped particle population quite 
differently from satellite absorption if the dust particles are small compared 
to the range of the protons (less than a fraction of a mm in our energy 
range). Strong forward scattering of light by both the E and G rings (Smith 
et al., 1982) suggests that these rings are largely made up of ym-sized 
particles. Under these circumstances, the ring absorption resembles energy 
loss in the upper atmosphere treated by Hess (1968), but the pitch angle 
dependence is quite different because atmospheric absorption occurs primarily 
near the mirror point and increases rapidly as the mirror altitude decreases. 
In contrast, at Saturn protons lose an energy aE each time they pass through 
the ring: 


AE 


w d^ 
cos 0 dx 


(19) 


where w is the mass per unit area of the ring (gm cm"2) , and dE/dx is the 
stopping power (MeV cm2 gm-1). This gives an average rate at which particles 
lose energy of (-dE/dt)- = 2AE/-cj3, if N(E)dE is the number of particles per 



21 


ORIGINAL PAGE IS 
OF POOR QUALITY 

In the Interval dE, then N(E) • (dE/dt) gives the flow of particles In 
energy, and the number lost Is given by the gradient. For a stable population, 
source Is required to cancel the loss (equ. 4.2 In Hess, 1968): 


4 §1 ■ s; (20) 

I 

where Sp 1s the source of particles. Solve equ. 20 and substitute the above 

I 

expression for dE/dt, j/V for N(E) and Sp/V for Sp, where V Is the proton 
velocity and Sp Is the proton flux Injected per second In units of cm“^ s”^ 
sr"^ MeV"^. The resulting expression for the differential proton spectrum Is: 


J 


Tp/r" cGse 
2w' TdrAOf)' 


f? 


/r 


dE' 


( 21 ) 


If \ni) are primarily Interested In the change of the energy spectrum due 
to absorption, then the energy dependence can be approximated by a power law: 
Sjj « E“^ , vp « and dE/dx « The value of 0.75 was calculated 

from the stopping power of silicon In the 10-160 MeV range and Is not much 
different for lighter elements In the rings. Substituting these values In 
equ. 21, we find 


j 


^ ^-(y-1.25) 


( 22 ) . 


The exponent y-1.25 shows that the proton spectrum should be significantly 
harder than the Injection spectrum and could turn over 1f y < 1.25. This Is 
In contrast to absorption by satellites where the proton spectrum should be 
softer than the local source. 
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Discussion 

The CRS observations measured the proton flux (Fig. 3) in two channels 
(48-63 and 63-100 MeV). Based on this data, we have determined the slope of 
the energy spectrum (Fig. 5) and obtained a general idea of changes in the 
angular distribution with dis*:ince (Fig. 6). Our analysis of the data 
demonstrated a substantial inward increase in the phase space density of 15 
BeV/G protons (Fig. 7) and thus confirmed the interpretation of Pioneer 11 
data by Van Allen et al. (1980). The most significant new result is the 
measurement of the energetic proton flux (63-160 MeV) in the Mimas absorption 
region. Because we can calculate the characteristic lifetime against 
absorption by Mimas, this measurement permits us to deduce the strength of the 
proton source required to maintain the observed flux. Based on the strength 
of the proton source, we can establish the strength of the GRAND neutron 
source., 

The particle distribution function in momemtum space obeys, on the time 
scale of interest to us, the steady-state diffusion equation (Schulz and 
Lanzerotti, 1974): 


3f _ 
~ 


,28 rO 3f 
8T 




(23) 


where D is the diffusion coefficient, 5 the local proton source and the 

p P ' 

loss term. Over the small range of L under consideration, we can consider D 
and D/L2 as constant and write equ. 23 asi 


D 


a^f 
— ? 



t = 0 
P 


(23a) 


If Tp is the characteristic lifetime then tp is given by: 
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Within statistical errors the distribution function is flat in the Mimas 

gap (Fig. 7); thus its second derivative vanishes or is at most small. 

Therefore by equ. 23a, the local loss and source terms balance in the gap. 

This intepretation is consistent with the discussion in the last section in 

which we showed that protons with pitch angles near 90° should be absorbed 

near the edge of the gap and that the flux in the gap has, therefore, to be 

replenished by a local source. In view of the non-relati vistic relation 

between the proton flux and distribution function, f = j/(2mE), and the short 

range in L, we can substitute the proton flux for the distribution function in 

equ. 23a and 24 to obtain Sp = j Ap where Sp was defined in equ. 7 and is 

equal to 2mE S . The required neutron flux is related to the proton source 

P 

strength by (Cooper and Simpson, 1980; B'^ake et al . , 1983): 


j 


n 



(25) 


where = lo^ sec is the neutron decay halflife, x the proton injection 
efficiency, and Jn the GRAND neutron flux. 

To calculate the integral proton source and neutron flux in the interval 
63-160 MeV, we calculated the effective energy with equ. 16c and used the 
value of to calculate Tp with equ. 15 (Table 3). Indications are that 
another small satellite or absorbing material exists in this region in 
addition to Mimas (Simpson et al., 1980; Vogt et al . , 1982), but this object 
is too faint to have shown up in Voyager images (Synott et al., 1981; Smith et 
al., 1982). Therefore we feel justified in using tp for Mimas because xp is 
proportional to r"2 a small object should not decrease xp signifi- 
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cantly. Values of x at L = 3.075 were obtained by interpolating between 
graphs of x vs. pitch angle calculated by Blake et al. (1983) for an isotropic 
neutron source. We measured a proton flux of 6 x 10-2 cm-2 s-1 sr-1 at a 45® 
Equatorial pitch angle inbound and 4 x 10“2 cm-2 s-1 sr-1 at 75® outbound. 

These and the other pertinent parameters are listed in Table 4. Based on these 
results, our best value for the neutron flux is 7 x 10-3 cm-2 s-1 sr-1. 

This value has been increased by 10%, as compared to the average of the two 
measurements in Table 4, in order to correct for protons above 160 MeV which 
were not counted. 

We have made a plausible argument that the source and loss terms in equ. 
23a should balance locally at the center of the Mimas gap; however, our 
counting statistics are too poor to prove rigorously that S^f/dL^ = 0. Should 
this not be the case, we have overestimated the required source strength. A 
lower limit for the source strength can be established by assuming that 
protons diffusing into the gap are uniformly distributed in the gap and that 
diffusion is loss free except in the Mimas absorption region. We integrate 
equ. 23a between 2.90 and 3.33 Rs* Because the first derivative of the 
distribution function vanishes at bo+h endpoints, we have: 

'-2 

/ § dL > / f/T dL (26) 

P Mimas ^ 

Over the small range in L, the proton source strength is approximately 
constant, and we ootain the following lower limit: 
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This proves that the required source strength can be no lower than 
30 percent of the value derived above (equ. 25). 

The experimental value of 7 x 10-3 cm-2 s-1 sr-1 for the neutron flux is 


consistent with calculations of the expected GRAND flux from Saturn's rings 
(Table 4). Cooper and Simpson (1980) estimated a neutron flux of - 10-2 cm"^ 
s-1 sr-1. Blake et al., (1983) obtained‘a value of 8 x 10-5 cm-2 s-2 for the 

source strength of protons above 10 MeV at 2 d65 R$ 3nd 0 «= 65°. Their value 
can be extrapolated to the orbit of Mimas by using the ratio of the injection 
efficiencies at the two distances. In terms of source strength per unit solid 
angle, the extrapolated value is 4 x 10-5 cm-2 s~2 sp~l as compared to our 
value of 9 X 10-7 cm-2 s-2 sr-1. Our value is significantly lower, but their 
figure also includes protons between 10 and 63 MeV. Saturn's atmosphere provides 
only a minor fraction of the GRAND neutrons because its flux is only < 10"^ 
cm-2 s-1 sr-1 (Cooper and Simpson, 1980), or more than 2 orders of magnitude 
below the required flux. 

Confirmation that the protons are produced by albedo neutron decay is 
offered by the proton spectrum and angular distribution. The slope of the 
spectrum above 48 MeV observed at 2.74 R5 resembles the GRAND proton spectrum 
at earth near L = 2.1 (Levine and Vette, 1970) and falls within the range of 
spectra predicted for a source at Saturn's rings (Blake et al . , 1983). A 
pancake distribution is also consistent with a neutron source at the ring 
plane provided the neutrons can escape isotropically. For this case, Blake et 
al. (1983) predict an angular distribution of approximately sinSe, which is 
close to our observations outside the absorption regions, that is between 3.5 R5 
and the orbit of Mimas and outside the G-ring to the orbit of Mimas (Fig. 6). 

Once the source strength is known, we can calculate the proton residence 
time in different parts of the inner magnetosphere. In a loss free diffusion 
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region this gives the time required to build up the observed flux, while it 
gives the lifetime against absorption in a region dominated by local losses. 
Table 5 lists both the omni-directional proton flux and the residence time 
which is given by the ratio between the flux and source strength 8.8 x 10“® 
protons cm-2 s"2j^ omni-directional flux at periapsis, 2.734 Rs» Is based 
on our measurement at a 55® pitch angle and sin^e distribution; ± 20 percent 
covers distributions from sin 0 to sin^e. The residence time of 30 y is 
consistent with those found in the inner radiation belt of the earth. The 
residence time is 12 y in the G-ring; this is the lifetime against absorption 
because both the first and second derivatives of the distribution function 
vanish in this region (Fig. 7b). The residence time has decreased to 1.4 y at 

3.33 Rg. Because the second derivative of the distribution function is 
negative, the local source is larger than the local loss term. We may 
conclude that diffusion away from 3.33 Rj plays a significant role in limiting 
the residence time. 

We are now able to compare omni-directional fluxes observed by different 
instruments. At 2.734 Rs» Pioneer 11 measured fluxes of 1.5 x 10^ cm-2 g-l 
(Van Allen et al , 1980) and 1.1 x 10^ cm-2 s-1 (Fillius et al , 1980) for the 

integral proton flux above 80 MeV. Using our energy spectrum (Ech ~ 

and the results given in Table 5, we obtain 5. 9 x 103 cm-2 s-1 above 80 MeV 

which is about half as much as the Pioneer 11 result. Krimigis (private 
communication) measured a flux of 3.8 x 10^ cm-2 s-1 of 16-160 MeV protons at 
the same distance but at a latitude of 8.2®; this translates to a flux of 4.5 
X 103 cm-2 s-1 at the equator for a sin^e distribution. This should be 
compared to our value of 8.2 x lO^ above 63 MeV. As can be seen, the 
different observations are within a factor of 2 of our value, and the spread 
can be accounted for by the substantial differences between the instruments 
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and the difficulty of obtaining a true omnl-dl recti o»jal flux In the presence 
of a very anisotropic pitch angle distribution. These observations are, 
therefore, consistent with the required residence time; however, they can not 
be used as proof of such stability. 

Local losses modify the observed proton spectrum as compared to the 
Injection spectrum. For losses due to absorption by dust rings or satellites, 
the direction of the change from the source spectrum can be deduced on the 
basis of the processes discussed In the previous section. The characteristic 
lifetime due to satellite absorption decreases with proton energy; therefore, 
spectra should be softer In the Mimas and Enceladus absorption regions. We 
were unable to measure spectral parameters at Mimas; however, a steepening of 
the spectrum was observed at Enceladus (Fig. 5). 

Rased on the model of dust-ring absorption, spectra should become harder In 
the fi and E-rlng. Specifically the Index y of a power law spectrum (j « E“^) 
should decrease by 1.? In the G-rIng where local absorption dominates. 
Unfortunately, the errors on our measuren»ent of y in this region are large, 
but relative errors should be substantially smaller, and the trend (Fig. 5) 
between 2.734 and 2.82 Rj is In the expected direction. This effect may also 
be the reason for the Inverse slope (negative y) observed by Krimigis and 
Armstrong (1982) between 20 and 50 MeV. The neutron spectrum should be almost 
flat In this energy range (Blake et al., 1983), and the slope of the proton 
spectrum should reverse. Any E-ring absorption is superimposed on effects due 
to Enceladus. Changes In the spectrum above 48 MeV Indicate that absorption 
by Enceladus dominates, while the 10-50 MeV data of Krimigis and Armstrong 
favor E-rlng absorption. It Is not yet clear how to reconcile these 
apparently contradictory results. 
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By now the reader may have become aware that we carefully avoided a 
discussion of the diffusion coefficient. Diffusion theory has been applied 
extensively to loss free diffusion and localized sources and sinks, but a more 
detailed analysis is required if extended source and loss terms are 
involved. Changes of the distribution function with radial distance have to 
be analyzed and the second derivative (equ. 23) is especially important. 

Such an analysis is planned to derive the proton diffusion coefficient and the 
density of the G ring. 

Conclusions 

The properties of the energetic proton flux (48-160 MeV) in the inner 
magnetosphere of Saturn are entirely consistent with a GRAND source at the 
rings. In particular, the residual flux in the sweeping region of Mimas is 
consistent with local injection from neutrons decaying in this region and 
requires a neutron flux of 7 x 10“3 cm-2 s-1 sr-1 above 63 MeV. This is 
consistent with the neutron flux expected from cosmic ray interaction with the 
rirgs (Blake et al., 1983; Cooper and Simpson, 1980). The angular distri- 
bution of the protons requires a nearly isotropic neutron source which 
suggests that the neutron absorption in the ring plane is relatively small. 

The proton trapping lifetime is comparable to that found in the earth's inner 
radiation belt. 

We can also ask whether our results can be explained equally well by 
episodic diffusion as the mechanism which transports protons past the 
satellites to populate the inner magnetosphere (McKibben and Simpson, 1980). 

The phase space density increases about 30 times between 4.3 and 2.74 ^ 

large diffusion mean free path during disturbed conditions can explain 
diffusion past Mimas but does not result in an increase in phase space density 
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as long as the diffusion process is reversible. This requires that the proton 
density outside of Mimas at 4.3 Rg must have been 30 times the observed value 
when the episodic diffusion mechanism operated. This problem does not arise, 
and the phase space density can increase inward if Inward displacements are 
always larger than outward steps (McKibben and Simpson, 1980). This would be 
the case if the mean free path is proportional to with m > 0; however, such 
irreversibility in magnetospheric diffusion has not yet been demonstrated. 

The proton spectrum provides further evidence; it is consistent with GRAND 
injection and considerably harder than would be expected from the inward 
diffusion of the proton population found in the outer magnetosphere. The 
variable and highly anisotropic angular distribution, proportional to sin^e 
with n = 4-7, would not be expected from inward diffusion. Any one of these 
points by Itself could be consistent with episodic diffusion, but collectively 
they strongly favor a GRAND source. 
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Appendl x 


Corrections for Counter Dead Time, Accidental Coincidences, 
and Particles Penetrating Through the Sides 


In the inner magnetosphere of Saturn the HET detector (Fig. 1) operated 
at much higher rates than normal, and extensive corrections were required to 
derive the true proton flux. These corrections depend on the time-constants 
and logics of the electronic system and on the new geometric configuration 
which was produced when detectors Bl, Cl, and the guard counters were 
deactivated. Three different corrections where required; (a) dead time 
correction for the counts lost, (b) subtraction of accidential coincidences 
to obtain true coincidence rates, and (c) corrections for higher energy 
particles that pass through the side of the detector. Because corrections to 
the B2 C4 TT rate were too large to obtain a valid flux between 27-48 MeV, we 
will discuss in the following only the B2 C4 C3C2 and B2 C4 C3 C2 coincidence 
rates. 

The logic of the HET coincidence-anticoincidence matrix is controlled by 
a strobe pulse. After the strobe pulse is triggered by a particle entering 
the telescope, the discriminators on each detector of the telescope are 
interrogated to see whether or not they were fired. The event is added to 
that rate counter which corresponds to the observed combination of triggered 
detectors; if no match is found, the event is discarded. The individual 
counting rates of each detector are monitored in space. Extensive ground 
calibrations permits us to derive dead time losses and accidential coincidente 
corrections from these rates. 


ORIGINAL PAGE IS 
OF POOR QUALnY 


31 


The dead time associated with each detector is controlled by the recovery 
time Tjj of the discriminator. If R is the true and N the observed rate of a 
detector, then 


TTVJTT 


(AU 


where is the fraction of time that the detector is sensitive. For HET 

Tjj ~ 16 MS with some dependence on pulse amplitude. The pulse height 
dependence was not included in this analysis, and we estimate that the 
uncertainty in x^j is ± 15%. Corrections calculated with equ. A1 were applied 
to coincidences involving the C2 and C3 detectors (Fig. 1). The combined dead 
time correction for the B2, C4, strobe (ST) coincidences was determined from 
large pulses which triggered the heavy ion slant threshold (SL). The trigger 
rate of SL was low enough that it required no correction. Most events cannot 
trigger SL without also triggering at least B2 and C4, and the ratio of the 
B2, C4, ST, SL rate to the SL rate gives the fraction of the large pulses 
which succeed in triggering B2, C4, ST coincidences. This ratio varied 
between 0.3 when B2 was 3 x 10^ counts s”^ and ~ 0.95 for 10^ counts s"^. 
However, it is not only a function of the B2 rate but also of the counting 
rates of other detectors in HET. Therefore, we attributed a ± 30% uncertainty 
to this correction. 

The calculation of accidental coincidences is simplified if we use the 
actual counting rates of the detectors rather than true rates because then the 
dead time insures that two counts cannot occur in rapid succession and be 
recorded as only one accidental count. Let x^ be the length of time that the 
strobe pulse will record a coincidence after a detector is triggered . For a 
counting rate N, and a true coincidence rate the fraction of the time that 
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accidental coincidences are possible is ta(N - Nt)‘* during this time strobe 
pulses not associated with true coincidences will give accidental coincidences. 
Therefore 


Ng = (ST - ST^) Tg (N - N^) (A2) 

where Ng is the accidentia! coincidence rate, ST the total strobe rate that 
can be in coincidence with N, and STt represents the strobe pulses in true 
coincidence. In our circuits tg depends on N and falls in the range of 26ys 
at low rates to 1 Bus at 1.5 x 10^ counts s-1. The accidental coincidence 
time was approximated by a quadratic in N and given a 15% uncertainty. 

Let us use the following notation for counting rates in applying equ. A2 
to NET: = B2 C4 'C3 . Ng = B2 CA C3 *C2 and N 3 = B2 C4 C3 C2 coincidence- 

anticoincidence combinations (the line over the symbol of a detector denotes 
anticoincidence). Use subscripts "a" and "t" to denote accidental and true 
coincidences and use the symbol of the counter for its counting rate. Then 
the strobe rate which can be in accidental coincidence with detector C2 is 
Nz + N 3 - N 3 f Therefore 


N33 = (Nj + Nj - N3t) T, CZ 


'^3 ' "3a ^ " 3 t 


(A 3 ) 


In the above expression we have ignored the small correction for true 
coincidences to the singles rate of detector C2. The equivalent expression 
giving the total accidental rate involving detector C3 is given by 
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N 


2a 


+ N 


3a 


[Nj + + Nj - + Nj^)] C3 


(A4) 


The resulting true coincidence-anticoincidence rates are therefore; 


No Tfl C2 

■^st ^3 1 - T, C2 

a 


N 


2t 


Ng Nj Tg C3 

1 - T, C2 1 - C3 

a u 


(A5) 

(A6) 


The dead time corrections discussed above were applied after the true 
coincidence rates had been calculated. 

Particle fluxes derived from these counting rates depend on the proton 
energy spectrum because the geometric factor is energy dependent. We did not 
use MET in one of its designed operating modes, and different path lengths 
through the stack were possible because the acceptance cone had been increased 
from 29° to 113°. Because the flux decreases with increasing energy, the 
upper energy cut-off is less important than the energy threshold which can be 
determined more accurately. The energy ranges and geometric factors shown in 
Table 1 are based on the following considerations. The lower energy threshold 
corresponds to just triggering the last detector of the stack for a particle 
trajectory at ^ 25° to the detector axis. The upper threshold corresponds to 
the highest energy proton with a high enough dE/dx to still trigger B2. 

In the quadruple coincidence mode (B2 C4 C3 C2) forward and backward 
moving particles fall into different energy ranges because of differences in 
shielding on the two ends. The geometric factor for forward moving particles 
is larger because these particles only have to reach the front surface of C2, 
while backward moving particles have to go through Cl or through the 
substantially thicker shielding of the side of C2. For triple coincidences. 
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differences in shielding account for the different energy ranges of particles 
that come in or go out through the sides. A Monte Carlo calculation of the 
geometric factor weighted by the particle spectrum would have given more 
accurate results but was not considered justified in view of the uncertainties 
associated with the dead time and accidental coincidence corrections. 

Fluxes in the energy ranges 48-63 and 63-160 MeV were calculated from 
triple and quadruple coincidences. Fluxes in the energy ranges given in Table 
1 were extrapolated from these two values on the basis of the assumed spectral 
shape. These fluxes were then used for calculating the effective geometric 
factor for quadruple coincidences and side corrections for triple 
coincidences. The calculations were performed reiteratively using power law 
(E-Y) or exponential (exp - E/Eqj^) spectra. The energy range covered is 
relatively small and the results of the two calculations agreed within a few 
percent. The errors associated with the 63-160 MeV flux were not increased 
sigificantly over those of the quadruple coincidence rate; however, because of 
uncertainties in the correction for high energy particles passing through the 
sides, the 48-63 MeV flux is considerably less accurate and untrustworthy near 
the Mimas absorption gap. 

As can be seen from Fig. 1, each of the coincidence conditions includes 
particles with substantially different path lengths through the detectors; 
therefore, the upper and lower 'energy thresholds are not sharp. The nominal 
48, 63, and 160-MeV threshold energies are not well determined and depend on 
the shape of the spectrum. Since few corrections have to be made to the 63 to 
160-MeV channel, the differential flux shown in Fig. 3 would change 
approximately proportional to the channel width, and the uncertainty in the 
width is expected to be less than 30 percent. Because corrections are 
required for particles penetrating through the sides, the effect of changing 
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the thresholds of the 48 to 63-MeV channel would be larger than for the 63 to 
160-MeV flux. A change in this flux would primarily affect the values of Eq^ 

and Y shown in Fig. 4. It is possible that the error bars shown beyond 3 R$ 
in Fig. 4 are optimistic; however, this does not affect the relative position 

of the points. A systematic error in Eqh would affect the relative phase 
space density shown in Fig. 7. The overall shape of the curve would not 
change but the magnitude of the increases and decreases could change. 
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Energy Range Geometric Factor 

(MeV) (cm2 sp) 


Quadrupl e/Coi nci dences 
of B2 C4 C3 C2 

Front incident 63 - 160 6.5 

Back incident 72 - 180 4,7 


Triple Coincidences 
of B2 C4 C3 with C2 anticoincidence 


Total geometric factor 

48 - 

63 

10.5 

Out through the side 


63 - 

160 

2 

In through the side 


72 - 

180 

2 

Detector Parameters: 

Detector 

Thickness 

(mm) 


Threshold Energy 
(MeV) 

Area 

(cm2) 

B2 

2 


2.16 

8 

C4 

6 


5.30 

9 

C3 

6 


5.28 

9 

C2 

6 


5.29 

9 
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Table 2 

Magnetic Cut-off Energies of Dipole Field 
(calculated using equations of Sauer, 1980) 


Dist. from 
Saturn (R^) 

Equatorial 

Field Line Distance* 

(Rs^ 

Latitude 

(Degree) 

Proton 
Energy (MeV) 

Alpha Particle 
Energy (MeV/n) 

in.o 

22.1 

180 

31 

7.5 


24.9 

300 

82 

5.0 

6.7 

28.9 

940 

302 

3.8 

5.0 

29.5 

1,940 

730 

2.74 

3.05 

20.0 

5,880 

2.570 

3.03 

3.05 

-9.9 

5,550 

2,410 

4.35 

5.0 

-19.5 

1,930 

708 

5.0 

6.1 

-22.6 

1,200 

402 

7.50 


-27.2 

265 

73 

10 


-28.6 

84 

22 


■*^From the "Z3" field model (Connerney et al., 1982) at distances where that 
model is applicable. These values correpond to the Equatorial distances 
shown in Figs. 2 to 7. 


ORIGINAL PAGE ii‘5 
OF POOR QUALITY 


42 - 


4/1 

n 

E 




o 

c 

cn 

to 

m 

«r>* 




E 

00 

cu 



• 

lO 

It 

CM 

■(S’ 

• 


• 


o 

o> 

rH 




J ^ 


CVI 

in 

CO 


o 

o 

in 


ft 00 


x: 

CO 

00 


o. 

’-3 

Ql 

r~“ 

X) 

ITS 


C 

o 

o 

L. 

a. 

>• 

0 ) 


o 


0 

E 

to 

40 



J 

J 

JZ 

JC 

ir> 

CO 

rH 

a> 

Q£ 

(X 



lO 

met 


» 

• 

LO 




•H 

• 

• 


o 


• 

in 

CO 


1—4 

CM 

CM 

II 

l-t 


o 


CM 

00 

r-H 







o 


m 

A 



CD 




« 

• 


CO 






CO 

o 






(O 

> 

fO 

<u 

•o 


<u 

u 

<u 

3 ; 


c 

o 


• • 

CL 

U 

0 

LT) 

B 

40 

40 

</> 

1 

J 

J 

J 

J 

x: 

CO 

o 

00 

CM 

00 

CT> 

00 




00 




to 


• 

• 

ir> 

o> 

o 

LO 

to 

lO 

• 

<11 


II 

O 

•«cf 

• 

o 

o 


o 


CTl 

CO 




rH 


o 

• 

in 


<n 

o 





<D 




o 

A 

A 

A 

A 



H- 

Im 

O 

•l- 






to 

OD 

rH 



irH 




l/> 

o 

u 

c 

fO 


TO 

c 

to 

4/1 

"O 

o 

(. 

CJ 

o. 


01 


K. 

o> 

CL 


o 


CO 


T3 

H 


c. 

(U 

4-’ 

o> 

E 

to 

i. 

a- 


i3 

H 


OJ 

> 


4-» 


<u 

c. 


TD 

O 


<u 

o. 

</) 

4 -> 

E 


/s 

<u 


Ol 


QC 


O) 




es» 

0 ) 


x» 

o 


<u 


o 

i_ 


p*. 


o. 


o 

/s 

a. 


0 ) 

E 


0 ) 


u 

&> 

> 

o 


E 

00 


C 7 I 

c 


w 


r:® 


M 


01 

c 

3 

0 
u 
c 
0) 

CSi 

L. 

<u 

01 

m 

& 

5 » 

<U 

xn 

•tJ 

O) 
c 
•»— 
L. 
(U 3 
I— TJ 

o> 

C to 
IQ 0) 
•U 

x: <o 
o c 

4-> -r- 
•r- -O 
Q. t- 
O 

f- O 
fO U 



c 





& 


•r* 


C 3 

(. in 


o 



Cft 


0) 


r“ 

■o 

o y 

O (O 


»f— 



c 


U 0) 



01 

44 S 

4-1 E 




4-> 

*»— 


c c 

40 

(A 

cn 

Ow 

<0 *r- 


Ob 


•»— 

Cl. 


<D O 

a 

U 

•D 

U 

3 £ 

•a 

u 


X3 

01 


4-> 

E 

•r* 

U (A 

CL CX 

cr 

o 

o 


U 

0) 


«/> C 

fm 

4-» I- 

01 0) 

<0 

Ui E 


40 


O 

3C 


*r- 

z 

Vi 

> r— 

<t- DTI 

o 

i. 

JQ 

X) 


tn 

CO 

T3 


•f— 

•0 o> 

O 

01 L. 

Of 

<C 40 

o 

4- 


3 

XI 

4#- 

&. un 

e 

01 

x: <«- 

a. 

40 

•r* 

O 

>*- 


»“ 01 

o 

4J 5*«» 

Si ^ 

01 je 

4J 


01 <ci 

l_ 


o 

TJ 

It? r“ 



B 

E 

(= 

<u 

DT> Ol 

01 

4/» 


<0 

U 01 

to 

U • 

«»“ x; 

•r* 

«A Ol 

o 

<C 

CL 

9 


C. 

•r* > 

3 

W 3 

u 

•M O 


c= 

u. i» 

o 

•f— 

+j 

o 

(Q 

•f" 

(. O' 

01 4J 

0) 4J 

« 

3 

01 01 

£_ 

•o 

TO 

L. 

X L. 

-o 

<0 01 

«*- -r- 

b- c 

*9 l~ 

O 

> a. 

>l 

CO 

•*— 

>> 

01 4J 

<0 

f ' 

•»- Q. 

«V» ifa 


CO 

«c 


ac 

tk 

C9 

>• 

QC 

O 


rri 

♦ <«■ 


0R1Q!IMAL . J 
OF POOR QUALiW 


T«blo 4 

Thfl GRAND Source Strength 
(63 to 160 MqV) end Theoretical Predlctlone 


Parameter Inbound Outbound 

Obaorved Proton Flux 6 x 10“^ 4 x 10"“^ 

(cm“^8“^8r“^) 

Local Pitch angle 67" 81" 

Equatorial Pitch angle 45* 75* 

X (injection afflclancy) 0.085 0.14 


Tp (characteristic llfotima) 

J/Tp (cm“^s“^Br“^) 

(Proton source strength) 

•^^n/XTp (cm"^8“^8r”^) 

(neutron flux) 

GRAND Proton Injftction at 65* 
above 10 MeV (cm“^e‘“^er'”^)* 

GRAND Flux from Rings 

(neutrons cm”^B~^Br“^)'^ 

GRAND Flux from Atmosphere 
(neutrons cm“*8~^ar'’^)’^ 


41 h 9 . 6 h 

4 X 10 "^ 1 X 10 “® 

5 X 10”3 8 X 10“3 

4 X 10 "^ 

~ 10~2 and < iQ-l 
10“^ 


Result of Rlake et al. (1983) extrapolated to the orbit of ^mae 
^ From Gooper and Slmpeon (1980) 
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Table 5: 

Omni-directional 63-160 MeV Proton Flux and Residence Time 


Equatorial Distance 

Proton Flux 

Residence Time 

(Rj) 

(cm-2s-l) 

(Approximate) 

2.734 

8.2 X 103* 

30 y 

2.80 (G Ring) 

3.2 X 103** 

12 y 

3.075 (Mimas) 


47 h+ 

3.33 

3.8 X 102++ 

1.4 y 


* Within ± 2.5 x 10^ for angular distribution from sine to sin^e 

** Calculated for the observed angular distribution of 0.4 + sine 
^ Actual lifetime against absorption calculated from absorption cross 
section of Mimas 

^ Calculated for a sin^e distribution 
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Figure Captions 

Figure 1 . Cross-sections of the B end of the High Energy Telescope > The SI 
detectors labeled A are 0.15 mm thick, B >■ 2 mm, and C 3 mm 
each. The Insulators were made of delrln (labeled DEL.), and 
screws that hold the telescope together are not shown. In Saturn's 
Inner magnetosphere, only detector B2 and the cross-hatched parts 
of the C detectors were active. Detectors Bl, Cl, and the guard 
rings around the C detectors were Inactive . 

Figure 2. Proton counting rates, 96s average, that have been corrected for 

accidental coincidences and loss due to counter dead time . The R 2 
rate, plotted at one-tenth of Its value, was due to B2, C4, C3 
coincidences with C2 in antl-colncldence; the R 3 rate was due to 
B2, C4, C3, C2 coincidences and accepted protons entering from 
either direction. The error bars shown at low counting rates 
represent statistical fluctuations, and at high rates they 
represent the uncertainties in the correction. The Z3 model of the 
magnetic field (Connerney et al., 1982) was used to determine the 
equatorial distance of field lines going through the spacecraft and 
to mark the expected position of satellite absorption features . 

The G ring, Mimas, and Enceladus are Identified by G, M, and Enc., 
respectively . 

Figure 3. The flux of energetic protons Is shown In two channels, 48-63 MeV 
and 63-160 MeV • The 48-63 MeV channel is not shown In the Mimas 
and G ring absorption region because of large uncertainties 
associated with corrections for protons entering or leaving through 
the side of detector C3. See caption to Figure 2 for a definition 
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of the distance scale and satellite Identification* The sharp dip 
near the G ring observed Inbound Just before peak flux was due to a 
spacecraft roll maneuver* 

Figure 4. The local pitch angles, solid line, were calculated relative to the 
observed magnetic field (courtesy of the Voyager Magnetometer 
Team) • Equatorial pitch angles shown as dashed lines were 
calculated from the local angles using the Z3 magnetic field model 
(Connemey et al* 1982). 

Figure 5. Indlcles of the proton energy spectra from the Inbound (solid dots) 
and outbound (open circles) passes ^ The top panel corresponds to a 
power law representation and the bottom panel to an exponential In 
energy. The equatorial field line distance Is based on the Z3 
model. 

Figure 6. The value of "n” Is shown versus equatorial field line distance for 
an assumed pitch angle distribution of the form sln^6. The ratio 
of the Inbound to outbound proton flux was used to calculate n; 
errors become very large near periapsls because of the small 
difference In the Equatorial pitch angle of the two observations. 

A 10-polnt running average was used between 4 .6 and 5 .5 Rg . 

Negative values correspond to a dumbbell distribution, zero to an 
Isotropic distribution, and positive values to a pancake 
distribution. Considering the size of the error bars, values 
beyond Enceladus are consistent with an Isotropic distribution. 

Figure 7 . Relative phase space density profiles for protons with constant 
first adiabatic invariant and 90® pitch angles. In view of the 
small change of the spectral Index with distance. It was kept 
constant with - 40 MeV for the 15 BeV/G curve and Eq^ “ 33 MeV 
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for the 60 BeV/G curve. Inbound values are plotted with crosses 

and triangles and outbound values with diamonds and circles. The 

2 

top panel was calculated assuming a sin 6 pitch angle distribution 
and the divergence between the Inbound and outbound curves shows 
that the angular distribution is changing. The bottom panel was 
calculated using a sln’^6 distribution where n changes with distance 
as shown in Figure 6, but local fluctuations had been smoothed 
out . This guarantees that the Inbound and outbound passes coincide 
and only the outbound pass has been plotted. Typical error bars 
are shown at different positions. It should be noted that the 
curves constitute only an upper limit beyond 5Rg because of the 
presence of a substantial cosmic ray component. 

Figure 8. Various particle absorption geometries are shown for protons with a 
gyroradlus p, a gyroperlod Tg, a drift velocity Vj, and a velocity 
component along the field line of Vj|. The radius of the satellite 
is denoted by r and the perpendicular distance between the path of 
the gyrocenter and the center of the satellite by x. (A) With 
X < p -■ r, the satellite moves inside the gyroclrcle and absorption 
occurs at the leading edge of the gyroclrcle only for gyrocenter 
locations along D; a symmetric situation occurs as the satellite 
emerges from the gyroclrcle. (B) With p-r < x < p + r, the 
satellite absorbs continuously between first and last contact over 
the distance 2D. (C) Protons that mirror off the equator and 

have the proper gyrophase can skip the satellite vertically if 


2r < V, T . 

I 8 








Fig. 3 
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